On the collapse in fourth order gravities 



B. Montes Nunez*, J. A. R. Cembranos 1 and A. de la Cruz-Dombriz** •* 

* Centra de Investigaciones Energeticas, Medioambientales y Tecnologicas (CIEMAT), 28040 

Madrid, Spain. 

' Departamento de Fi'sica Tedrica I, Universidad Complutense de Madrid, E-28040 Madrid, Spain. 
** Astrophysics, Cosmology and Gravity Centre (ACGC), University of Cape Town, Rondebosch, 

7701, South Africa. 

* Department of Mathematics and Applied Mathematics, University of Cape Town, 7701 
Rondebosch, Cape Town, South Africa. 

Abstract. The gravitational collapse in fourth order theories of gravity defined by an arbitrary 
action of the scalar curvature shows significant deviations with General Relativity. The presence of 
a new scalar mode produces a higher initial contraction that favors the reduction of the collapsing 
time. However, depending on the particular model, there are fundamental differences when the 
modifications to the General Relativity collapse leave the linear regime. These analyses can be 
used to exclude an important region of the parameter space associated with alternative gravitational 
models. 
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INTRODUCTION 

The gravitational collapse for a spherically symmetric stellar object has been analyzed in 
detail in the standard General Relativity (GR) theory [1]. By assuming the metric of the 
space-time to be spherically symmetric and that the collapsing fluid is pressureless, the 
found metric interior to the object turns to be Robertson-Walker type with a parameter 
playing the role of spatial curvature and proportional to initial density. The time lapse 
and the size of the object are given by a cycloid parametric equation with an angle 
parameter y/. 

In spite of the fact that GR has been one of the most successful theories of the 
twentieth century, it does not give a satisfactory explanation to some of the latest 
cosmological and astrophysical observations with usual matter sources. On the one hand, 
the baryonic matter content has to be supplemented by a dark matter (DM) component. 
Although there are many possible origins for this additional component [2], DM is 
usually assumed to be in the form of thermal relics that naturally freeze-out with the right 
abundance in many extensions of the standard model of particles [3]. Future experiments 
will be able to discriminate among the large number of candidates and models, such as 
direct and indirect detection designed explicitly for their search [4], or even at high 
energy colliders, where they could be produced [5]. 

On the other hand, a dark energy contribution needs to be considered to provide 
cosmological acceleration [6]. Alternatively, a modified gravitational interaction can 
potentially generate the present expansion of the Universe without invoking the presence 
of any exotic dark energy [7]. In this context, fourth order theories of gravity defined by 



adding to the action a function of the scalar curvature: 



has been extensively studied during the last years [8, 9, 10, 11]. 



GRAVITATIONAL COLLAPSE 

In order to study the gravitational collapse inside this modified gravity theories, so-called 
f(R) theories, we can make use of the spherically symmetric metric. The collapsing 
object is considered to be approximated pressureless and homogeneous. Therefore, we 
can search a separable solution for this metric as follows [1]: 

ds 2 = dt 2 - U(r, t)dr 2 - V (r, t) (d9 2 + sm 2 9d(j) 2 ) , (2) 
U(r,t)=A 2 (t)h(r), V(r,t)=A 2 (t)r 2 . (3) 

Components tt, rr and 99 for the modified tensorial equations may be written re- 
spectively in terms of the functions A{t) and h(r). Once we have calculated h(r), 
the resulting metric is formally the same as the one obtained in GR. We consider 
p(t) = p(t = 0)/A(t) 3 (given by the energy motion equation for dust matter). Fur- 
thermore, provided that the fluid is assumed to be at rest for t = 0, initial conditions 
A(t = 0) = and A(t = 0) = 1 hold. We define R(t = 0) = R and p (t = 0) = p to 
simplify notation. Taking into account the previous conditions, the evolution equation 
for A is found to be [12]: 

A2 = ~ Zn Ttt ,A ^Gpo(2-MRo))-f(Ro)(l+f R (Ro)) (4) 

-3f R (Ro)f R (Ro)} + (7^72) ^pPo^ 1 ~ 6(1 1 /2) [*xG Po A- l f R 

+3A 2 f R f R -3AAf R (l-f R ) +A 2 f(R)(\ +f R )] . 

Let us remind at this stage that the zeroth order solution of GR is given by the parametric 
equations of a cycloid [1]: 

w + sin \if 1 . 

t= 2Vk A G = -(l+cosv). (5) 

Expression (5) clearly implies that a sphere with initial density po and negligible pres- 
sure will collapse from rest to a state of infinite proper energy density in a finite time 
that we will denote Tq. This time is obtained for the first value of y such as A<j = 0, 
i.e. for y = n. In order to study the modification to the gravitational collapse in f(R) 
theories, we have expanded A around Aq and f(R) around the scalar curvature in GR 
(R = R G ) [12],i.e. A = A G +g(\if), and f(R) ~ f(R G ) + f'(R G )(R-R G ). The presence of 
a function f(R) in the gravitational Lagrangian represents a correction of first order with 
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Figure 1. Validity of the perturbative regime for f(R) = sR 2 [8, 1 1] (left panel) and f(R) = eR 1 [9] 
(right panel). As it is indicated in the figures, in the first case, it is the high density region, the one that 
exhibits a non linear regime; whereas it is the low density region for the second model (on the right). The 
right area on the left panel is excluded by the condition e < 2.3 x 10 22 GeV -2 [11]. Finally, the density 
marking the beginning of structure formation (SF) Psf — 1.5 x 10~ 38 GeV 4 , and the dark energy (DE) 
density p^E — 2.8 x 10~ 47 GeV 4 , have also been plotted for reference [12]. 



respect to the usual Einstein-Hilbert action. By substituting the above series expansions 
in expression (5) until first order in £, we find that equation (5) becomes [12]: 

ig(f)s' = ~oos- 2 (^)g + ^cos 2 {^'j(f(Ro ll )+3kfR(Rco)) (6) 

- >°> - Ik™* (!) tW> +^f in (I) cos* (|) M *o) ■ 

By studying this equation for particular f(R) models, it is possible to constraint an 
important region of the parameter space of this gravitational theories (See Fig. 1). This 
validity tests are interesting since the higher contraction that these models generally 
present in the perturvative regime, can alleviate some tension of the standard ACDM 
model with observations of structures at high redshift [13]. 
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